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^ ' Recent observationsBi of high-redshift supernovae seem to suggest that the 

^ ■ global geometry of the Universe may be affected by a 'cosmological constant', 
t:;^- ■ which acts to accelerate the expansion rate with time. But these data by them- 
O , selves still permit an open universe of low mass density and no cosmological 
Q"^ ! constant. Here we derive an independent constraint on the lower bound to the 
^ \ mass density, based on deviations of galaxy velocities from a smooth universal 
expansiontm. This constraint rules out a low-density open universe with a van- 
O \ ishing cosmological constant, and together the two favour a nearly fiat universe 
' in which the contributions from mass density and the cosmological constant are 

■ comparable. This type of universe, however, seems to require a degree of fine 
^ ■ tuning of the initial conditions that is in apparent conflict with 'common wis- 

^ . dom'. 

The standard cosmological model based on Einstein's gravity and global homogeneity is 
characterized by two fundamental parameters that measure the main contributions to the 
total energy density (in terms of a 'critical' density): the mean mass density and the 
cosmological constant i^A. The latter, an intrinsic part of the theory of general relativity, 
represents a uniform energy density that is associated with the vacuum (not the mass), and, 
if positive, acts as a repulsive component of the global gravitational force field. The values of 
these parameters determine the type of the universe we live in, its global extent and ultimate 
fate, as follows. The parameter plane (Fig. 1) is divided by three lines into six regions. The 
line fijn + = 1 defines a 'flat' euclidean geometry, favoured by theories of inflation in 
the early universe; it separates curved-space (non-euclidean) models of 'closed' finite extent 
(above) and 'open' infinite extent (below). The horizontal line f^A = roughly distinguishes 
between an 'unbound' universe that will expand forever (above) and a 'bound' universe 
that will eventually re-collapse (below). The special point (fimj^A) = (1)0) corresponds to 
the simplest model termed Einstein-deSitter; unless the universe exactly fits this model, it 
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evolves away from this point as it expands. If ^Ia > 0, the universe eventually crosses the line 
Qa = 0.5Qrn, which separates between a universe whose expansion decelerates by attraction 
(below) or accelerates by repulsion (above). 

New techniques to measure distances of supernovae (type la) at large distances (cosmo- 
logical redshifts 0.5 — 1), where the large-scale curvature of space-time plays a noticeable 
role, enable a classical cosmological test based on how the time dependence of the Hubble 
relation between velocity and distance depends on the cosmological parameters!. The corre- 
sponding, inclined (blue) confidence limits in Fig. 1 are based on the results of the Supernova 
Cosmology Projectu, which are fully consistent with the findings of the High-z Supernova 
Search Teamll. The allowed region in the parameter plane is an elongated stripe, crudely 
approximated by 0.8f2m — O.Gf^A = — 0.2±0.1 (la errors). Based on this result alone, a nearly 
flat universe is likely, with a broad maximum to the probability near (f2m,f2A) = (0.3,0.7), 
but an open model, of (0.1,0), for example, is still allowed at the 2% confidence level. An 
orthogonal constraint is required in order to remove the degeneracy between Qa and fim- 

Here we show that such a constraint is provided by current data of 'peculiar' velocities 
of galaxies on scales of ~300 million light years. These velocities, which are produced by 
gravity due to local fluctuations about the mean mass density, depend also on the value of the 
mean density itself, Q^. The constraints on from these velocities are almost independent 
of ^A, but combined with the inclined supernovae constraints any bound on fim effectively 
serves as a bound on ^a- 

We use two recent catalogues of galaxy peculiar velocities. One, named Mark con- 
tains distances for ~3,000 galaxies within ~70h~^Mpc; the other, named consists 
of ^ 1,300 spiral galaxies with a more uniform sampling in a similar volume. In order to 
obtain the pecuhar velocity of a galaxy, its total velocity is measured using the Doppler 
redshift, and its distance is separately inferred using the so-called TuUy-Fisher method, with 
an error of 15 — 21%. The desired peculiar velocity along the line of sight is obtained by 
subtracting, from the total velocity, the Hubble expansion velocity, at the galaxy distance. 
The peculiar velocities are carefully corrected for systematic errorsEHlHl. 

The POTENT methodlll'lii assumes that the peculiar velocities were generated by gravity 
and that they therefore represent a potential flow on large scales. This allows a recovery 
of the three-dimensional velocity field of galaxies in our local cosmological neighbourhood. 
For an assumed value of fl^, the underlying field of mass-density fluctuations can be ex- 
tracted, using slightly nonlinear approximations to the relation between velocity and density 
fluctuations within the theory of gravitational instability. Assuming also that the initial 
fluctuations were a gaussian random field, direct dynamical constraints on can be ob- 
tained without appealing to the observed spatial distribution of galaxies. These constraints 
are thus independent of the unknown relation between galaxies and the underlying mass 
density (termed 'biasing'). Several different methods applied to the POTENT Mark III 
fields consistently yield a lower bound of Qrn > 0.3 at the 99% confidence level, marked by a 
vertical red line in Fig. 1. One method constrains fim based on the large diverging outflow 
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Figure 1: Constraints in the 0,^ — J^a plane, showing relative confidence limits of 68%, 90% and 
99%. The inclined contours (blue, same in both panels), that roughly constrain O.Srim —Jl).6Q\, 
arise from the global geometry of space-time based on supernovae as distance indicatorall. The 
lower bound > 0.3 (red) represents a 99% confidence limit based on several different studies 
of peculiar velocities, using POTENT reconstruction, independent of the biasing relation between 
galaxies and masaJ^I. The almost vertical contours (red) arise from the peculiar velocities of the 
Mark III (left) and SFI (right) catalogues, based on a likelihood analysis which assurojes a parametric 
mass power spectrum of the CDM family with COBE normalization on large scaleda'Q. The central 
contours shown are for a standard case of fixed h = 0.65 and n = 1; the uncertainties in these 
values are represented by the outer 99% contours on both sides (long dashed red), computed for 
the extreme cases of fixed (/i, n) = (0.55,0.9) and (0.75,1.1). These encompass a conservative 
range of non-negligible likelihood based on the flow data. The strong bounds from flows are on 
r^ni; the apparent upper bound on f^A arises indirectly from the COBE normalization and is very 
uncertain; the contour can extend further in the Oa direction, allowing significant overlap with the 
supernova confidence region. The corresponding joint (relative) confidence limits from supernovae 
and flows are shown (black contours); they are almost the same for Mark III and SFI. The combined 
constraints thus favour an unbound and roughly flat universe with comparable contributions from 
the cosmological constant and mass density. A universe with very low and zero cosmological 
constant is ruled out. 
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recovered in the vicinity of a large nearby underdense region (void) , using the fact that even 
an empty void cannot induce outflows as large as those observed if the density outside it is 
too IoweI. a second method uses the deviations of the initial fluctuations from a gaussian 
distribution, when derived from the velocity data, assuming a value of f2m that is too low§. 
A third method derives Q^a based on the gravit at ionally- induced deviations from gaussian 
distribution of the present-day velocity-divergence fielJ. Similar constraints confirming the 
lower bound of f2m > 0.3 are currently being obtained from the SFI data. 

In a new analysisi'0, not based on POTENT, we have applied a maximum-likelihood 
analysis to the Mark III and SFI peculiar velocities to determine cosmological parameters. 
It is done via a parametric model for the mass-density fluctuation power spectrum P^, which 
measures the distribution of power among different length scales, proportional to k~^. This 
analysis uses linear gravitational theory and assumes that both the fluctuations and the 
errors are gaussian. The model used is of the popular cold dark matter (CDM) model 
family; is normalizedlli by the fluctuations in the cosmic microwave background (CMB) 
as measured on very large scales by the COBE satellite. The result provides constraints on 
permissible combinations of three parameters: Q^n, the power index n (where Pk (x on 
large scales), and the Hubble expansion constant h. These parameters enter via the shape 
of Pk as well as from the geometry and dynamics of space-time. The constraints obtained 
from the two data sets turn out to be similar; they define a two-dimensional surface in the 
Q^-h-n space. In the case of a flat cosmology this surface can be crudely approximated by 

^65^ = 0.58 ± 0.12 (where /igs is the Hubble constant in units of 65kms~^Mpc~^, and 
the errors refer to 90% uncertainty). 

The obtained Pk is determined by the velocity data (on scales of ~50 h~^Mpc) and is not 
much affected by the COBE normalization (at ~1, 000h~^Mpc); a similar Pk is reproduced 
to within a few percent when the amplitude on large scales is left as a free parameter. Also, 
the result is not dominated by the assumed distance errors; Pk is reproduced to better than 
10% when an error model with free parameters is incorporated in the likelihood analysisi. 
This robustness is within the la limits of the likelihood analysis. 

The peculiar velocities have also been analysed by methods that incorporate the spatial 
distribution of galaxies. For example, a comparison of the Mark HI velocities and the IRAS 
1.2Jy redshift surveylil yields results roughly consistent with the CDM model of fim ~ 0.5 
favoured by the Pk analysis of flows reported above. Our current bounds on from flows 
also partly overlap with constraints from the evolution of galaxy cluster Jll (fim = 0.45 ±0.2). 
However, a comprehensive joint analysis of all the available constraints is beyond the scope of 
this work. In particular, a thorough interpretation of the results involving redshift surveys 
(which span a noticeable range) must include a full discussion of the non-trivial biasing 
relation between galaxies and mass0, and of the complex systematic effects involved in these 
different studies. Instead, we choose to focus here on one set of constraints, the dynamical 
constraints from flows, combined with the supernova constraints. 

In order to obtain the desired constraints in the i^m-^A plane, we have applied the Pk 
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likelihood analysis to the Mark III and SFI velocities, using the COBE-normalizedEl CDM 
models with varying fim and that now span the whole parameter plane. The velocities 
are expected to be insensitive to the value of I^a^, but a weak dependence occurs as a 
result of the COBE normalization imposed. This is responsible for the slight left bending of 
the red vertical likelihood ridges (Fig. 1), and the apparent upper bounds on Q\. The la 
error in the COBE normalization translates to additional uncertainties of ±6% in and 
±20% in ^A, which are on the order of the uncertainties displayed by the likelihood contours 
shown. Thus, the COBE errors would tend to stretch the likelihood contours vertically along 
the Qa axis, weakening the apparent upper bounds on ^a- 

Since the velocity data favour an extended two-dimensional surface in the Qrn-h-n space, 
the likelihood analysis cannot determine the three parameters simultaneously; two of them 
should be fixed initially. We adopt as our standard case the simplest, scale-invariant {n = 1) 
initial spectrum, and h = 0.65 (as favoured by nearby supernova dataS). We crudely 
estimate a ±15% uncertainty in the Hubble constantlll, and ±10% in the power indexil. 
In order to illustrate the sensitivity of the constraints to these parameters, we have also 
determined fl^ and Q\ for two extreme cases of h and n values. The outer 99% contours 
are shifted accordingly in Fig. 1. 

The likelihood ridge of SFI extends into higher values of Q\ than for Mark III. The 
independent constraints from flows (using Pk) and supernovae overlap significantly for SFI, 
but only near the 90% contours for Mark III with {h,n) = (0.65, 1) (which improves when 
COBE errors are considered or if ri > 1 or /i > 0.65). The large uncertainty in the upper 
bound on ^a indicates that the degree of overlap should not be interpreted as indicating 
either consistency or inconsistency between the velocity and supernova data and the assumed 
standard cosmological model. 

A joint parameter estimation from supernovae and flows is therefore meaningful, although 
it is limited to relative likelihoods; a measure of absolute probabilities is not straightforward, 
despite a seemingly acceptable goodness of fit. The black joint contours (Fig. 1) are computed 
by multiplying the two likelihood values at each point, having assumed that the data are 
independent. The most likely joint values (shown as crosses in Fig. 1) for the supernovae 
and the combined Mark III and SFI dataset are {flm, ^a) = (0.5, 0.8), while the conservative 
joint 99% confidence limits allow values in the range 0.3 — 0.9 and Qa values in the range 
0.1 — 1.4. The constraints obtained separately from Mark III and SFI are very similar. 

The constraints from local flows thus remove the degeneracy in the constraints from the 
global-geometry test based on supernovae (and vice versa), and help rule out an open model 
with zero cosmological constant. A nearly fiat universe, with comparable contributions of 
matter and cosmological constant to the total energy density, is likely. The favoured model 
is therefore an unbound universe that will accelerate forever, although it is impossible yet 
to determine whether the global geometry is fiat, open or closed. 

Such comparable contributions from the mass density and the cosmological constant 
represent a puzzling fine tuning, for example, because the two parameters, Q^n and ^a, are 
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expected to vary with time in opposite senses. The standard theory expects the cosmological 
constant either to vanish or be larger by many orders of magnitudeiil. Perhaps entropic 
arguments may be needed to explain such a fine tuning. Although the observed constraints 
are not yet finally confirmed, they already seem interesting enough to pose a serious challenge 
to theoretical physics. 

Other constraints in the fim-^A plane are worth mentioning in this framework. Con- 
straints consistent with the supernova ridge but of larger uncertainty arise from the age of 
old star clusterS versus the Hubble expansion rate. Constraints of orthogonal orientation, 
roughly on l^A + ^m, can be deduced from the acoustic peaks in the sub-degree angularpower 
spectrum of fiuctuations in the CMB, as observed from balloons and from the groundiliHi. 
Two CMB satellites planned for the next decade, MAP and Planck, are expected to provide 
more accurate constraintsil'0. 

Future peculiar-velocity data are expected to improve the accuracy of the constraints. In 
addition to many new velocities based on distance indicators of the Tully-Fisher type, the 
most promising sources of large-scale peculiar velocities in the future will probably be the 
siblings of the same supernovae discussed above, but at low redshifts; their distances can be 
measured with 5 — 10% accuracy out to large distances and their sampling density is limited 
only by the patience of the observers. 
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